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Optimization  over  multi-order  cones 

Baha  M.  Alzalg*  and  K.  A.  Ariyawansa^ 


Abstract 

In  this  paper  we  propose  multi-order  cone  programs  (MOCPs)  as  a  new  class  of  convex 
nonlinear  optimization  problems  that  includes  linear  programs,  (convex)  quadratic  programs, 
second-order  cone  programs  and,  more  generally,  p  -order  cone  programs  as  special  cases.  In 
MOCPs  we  minimize  a  linear  objective  function  over  the  intersection  of  an  affine  set  and  a 
product  of  multi-order  cones.  We  refer  to  them  as  deterministic  multi-order  cone  programs 
(DMCOPs)  since  data  defining  them  are  deterministic.  We  present  the  definition  of  DMOCPs 
in  primal  and  dual  standard  forms.  Then  we  introduce  two-stage  stochastic  multi-order  cone 
programs  (SMOCPs)  (with  recourse)  to  handle  uncertainty  in  data  defining  DMOCPs  and  de¬ 
terministic  mixed  integer  multi-order  cone  programs  (DMIMOCPs)  to  handle  DMOCPs  with 
integer-valued  variables.  We  describe  an  applicational  setting  and  present  DMOCP,  SMOCP, 
and  DMIMOCP  models  arising  in  that  setting. 

Keywords:  Linear  programming;  Stochastic  programming;  Recourse;  Second-order  cone  pro¬ 
gramming;  Mixed  integer  programming 


1  Introduction 

Semidefinite  programming  [16,  13]  problems  were  extensively  studied  during  the  late  1990s  as  a 
class  of  optimization  problems.  They  are  extensions  of  linear  programs  and  provide  novel  modeling 
capabilities.  Interior  point  algorithms  could  be  derived  for  them  (often  utilizing  their  symbolic 
similarities  to  linear  programs). 

Ariyawansa  and  Zhu  [5]  (see  also  [10])  presented  stochastic  semidefinite  programs  that  extended 
stochastic  linear  programs  [17,  6],  and  allowed  the  derivation  of  elegant  interior  point  algorithms 

?!•  1°J- 

It  soon  became  apparent  [9,  1]  that  almost  all  applications  of  semidefinite  programs  indeed 
lead  to  a  subset  of  semidefinite  programs  termed  second  order  cone  programs.  We  refer  to  them  as 
deterministic  second  order  cone  programs  (DSOCPs)  because  they  are  defined  using  deterministic 
data.  In  DSOCP  we  minimize  a  linear  function  over  the  intersection  of  an  affine  set  and  a  Cartesian 
product  of  second  order  cones. 

In  this  paper,  we  present  three  extensions  of  DSOCPs.  First,  we  present  primal  and  dual  forms 
of  (deterministic)  multi-order  cone  programs  (DMOCPs)  in  which  we  minimize  a  linear  function 
over  a  Cartesian  product  of  pth- order  cones  (we  allow  different  p  values  for  different  cones  in  the 
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product).  We  present  generic  applications  that  extend  those  in  [1,  Section  2.2],  We  also  present  a 
glimpse  of  the  duality  theory  that  we  are  developing  in  [2]  for  DMOCPs. 

Second,  we  present  two-stage  stochastic  multi-order  cone  programs  (SMOCPs).  SMOCPs  are 
a  way  of  handling  uncertainty  in  data  defining  DMOCPs.  Then  we  demonstrate  that  stochastic 
linear  programs  and  stochastic  quadratic  programs  are  special  cases  of  SMOCPs. 

Our  third  extension  is  introduced  to  handle  modeling  situation  in  which  some  of  the  variables  in 
an  optimization  problem  can  only  take  integer  values,  or  even  0  or  1.  This  leads  to  (deterministic) 
mixed  integer  multi  order  cone  Programs  (DMIMOCPs)  and  0-1  deterministic  multi  order  cone 
programs  (O-lDMOCPs). 

We  then  demonstrate  how  decision  making  problems  associated  with  facility  location  problems 
lead  to  a  DMOCP  model,  an  SMOCP  model,  a  O-IDMOCP  model,  and  a  DMIMOP  model. 

We  begin  with  an  introduction  to  our  notation. 


1.1  Notations 

We  begin  by  introducing  some  notations  we  use  in  the  sequel. 

Let  Rmxn  and  RnVn  denote  the  vector  spaces  of  real  m  x  n  matrices  and  real  symmetric  n  x  n 
matrices  respectively.  For  U,  V  £  RnVn,  we  write  U  P  0  (U  y  0)  to  mean  that  U  is  positive 
semidefinite  (positive  definite),  and  U  >z  V  or  V  A  U  to  mean  that  U  —  V  y  0. 

All  vectors  we  use  are  column  vectors  with  superscript  T  indicating  transposition.  We  use 
for  adjoining  vectors  and  matrices  in  a  row,  and  use  for  adjoining  them  in  a  column.  So,  for 
example,  if  x,  y,  and  z  are  vectors,  we  have: 


x 


y  |  =  (*  ,V  ,z  )  ={x;y;z). 
z 

If  A  C  and  B  C  Rz,  then  the  Cartesian  product  of  A  x  B  :=  {(cc;  y)  :  x  £  A  and  y  £  £>}. 
For  each  vector  x  £  Rfc  indexed  from  0,  we  write  x  for  the  sub-vector  consisting  of  entries  1 
through  k  —  1;  therefore  x  =  (xq;  x). 


Given  p  >  1,  the  j//l-order  cone  of  dimension  n  is  defined  as  Q”  :=  {x  =  (xq',x)  £  R  x 


pn— 1 


xq  >  | \x\ |p}  where 


denotes  the  p-norm.  The  cone  Q1  is  convex,  pointed,  closed  and  with  a 


nonempty  interior  (see,  for  example,  [18]).  As  special  cases,  when  p  =  2  we  obtain  the  second-order 
cone  (also  known  as  the  quadratic,  Lorentz,  or  the  ice-cream  cone)  of  dimension  n,  and  when  p  =  1 
or  oo,  Qp  is  a  polyhedral  cone. 


We  write  x  A  0  to  mean  that  x  £  Q” ,  and  x  P  y  to  mean  that  x 


Given  1  <  p*  <  oo  for  i  =  1,2,' 


,  r.  Let  Q 


write  x  0  to  mean  that  x  £  Q^1’”2  , 

-{Pl,P2,---,Pr}  \P1,P2,—  iPr) 

.  (ni,n2,-,nr)  n 

y  —{pi,P2,  "  ,Pr) 

It  is  immediately  seen  that,  for  every  vector  x  £  R’ 

and  only  if  x  is  partitioned  conformally  as  *  =  (cci;  X2 ; 1 
For  simplicity,  we  write: 

M) 


(ni,n2,-"  ,nr) 
(P1.P2,-  ,Pr) 


y  tg  "• 


a. 


Pr  * 


:=  Qni  X  Qn2  x  •• 

^P1  ^P2 

’71r'>  and  x  y  to  mean  that  x 

—  \Pl,P2,---,Pr)  ^ 


We 


where  n  =  Yli=i  no  x  (£ 
_{m) 


(ni,n2,-"  ,nr) 
{Pl,P2,-  ,Pr) 


0  if 


;  xr)  and  aq  0  for  i  =  1, 2,  •  •  •  ,  r. 


•  Qp  as  Qp  and  x  [  (1  as  x  (Z{p)  0  when  n  is  known  from  the  context; 
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Qti*2-'-Z)  as  Q<PuP2,-,Pr)  and  *  0  as  *  h{pi:P2,...,Pr)  o  When  m,n2,' 

are  known  from  the  context; 


_(ni,n2,'"  !»r) 


•  x  >~ 


,P) 


0  as  x  y(n]  0; 


x  y  0  as  x  y  0;  and  x  <Zr(2)  0  as  x  yr  0  when  the  problem  includes  only  (linear  and) 
second-order  cone  constraints. 


Note  that,  for  every  vector  *  G  Mn  =  Rfc  x  M.k  x  •  •  •  x  R  ,  x  y.r^  0  implies  that  x  is  partitioned 
regularly  as  x  =  (, x±,X2 ; \xr )  and  each  subvector  Xi  lies  in  the  pth- order  cone  of  dimension 
k  =  n/r  for  i  =  1, 2,  •  •  •  ,  r.  In  this  case,  if  n  =  r  and  p  =  2,  then  x%  €  Q\  =  {t  €  M  :  t  >  0}  for  each 
i  =  1,  2,  •  •  •  ,n.  So  x  y,n  0  means  the  same  as  x  >  0,  i.e.,  x  lies  in  the  nonnegative  orthant  of  R". 


2  Definition  of  a  DMOCP 

Let  r  >  1  be  an  integer,  and  Pi,P2,-"  ,Pr  are  such  that  1  <  Pi  <  oo  for  i  =  1,2,---  , r*.  Let 
m,  n,  ni,  n2,  ■  •  •  ,  nr  be  positive  integers  such  that  n  =  Then  we  dehne  a  DMOCP  in 

primal  standard  form  as 

min  cTx 

( P )  s.t.  Ax  =  b 

x  y ,  ,0 

—  (pi,P2,-,Pr) 

where  A  €  Rmxn,  b  G  Rm  and  c  G  Rn  constitute  given  data,  *  G  Rn  =  Rni  x  Rn2  x  •  •  •  x  Rnr  is 
the  decision  variable.  We  define  a  DMOCP  in  dual  standard  form  as 

max  bTy 

(D)  s.t.  AJy  +  z  =  c 

z  y,  ,0 

—  (91,92,  -,qr) 

where  y  G  Rm  and  z  G  Rre  =  Rni  X  Rn2  x  •  •  •  x  R"r  are  the  decision  variables,  qi,q2,  ■  •  •  ,qr  are 
integers  such  that  1  <  qi  <  oo  for  i  =  1,  2,  •  •  •  ,  r. 

If  (P)  and  (D)  are  defined  by  the  same  data,  and  qt  is  conjugate  to  pt,  in  the  sense  that 
1  jpi  +  1/qi  =  1  for  i  =  1, 2,  •  •  •  ,  r,  then  we  can  prove  relations  between  (P)  and  (D)  (see  Subsection 
2.2)  justify  referring  to  (D)  as  the  dual  of  (P)  and  vice  versa. 

2.1  Special  cases  of  DMOCPs 

A  deterministic  pth-order  cone  programming  (see  also  [12])  (DPOCP)  problem  in  primal  standard 
form  is 

min  cTx 

s.t.  A  x  =  b  (1) 

x  y  /  \  o 

—rip) 
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where  m,  n,  n\,  n,2,  ■  ■  ■  ,nr  are  positive  integers  such  that  n  =  E[=in*>  P  ^  [l,°o],  A  G  Wnxn, 
b  G  and  c  G  Rn  constitute  given  data,  x  G  Rn  is  the  primal  variable.  Clearly,  DSOCPs  are  a 
special  case  of  DMOCPs  with  Pi  =  p  >  1  for  all*  =  1, 2,  •  •  •  ,  r.  According  to  (D),  the  dual  problem 
associated  with  DPOCP  (1)  is 

max  bTy 

s.t.  AJy  +  z  =  c  (2) 

z  y  , >  o 

r{q) 

where  y  G  Rm  and  z  G  Rn  are  the  dual  variables  and  q  is  conjugate  to  p. 

Deterministic  second-order  cone  programs  (DSOCPs)  are  a  special  case  of  DPOCPs  (and  hence 
of  DMOCPs)  which  occurs  when  p  =  2  in  (1)  (and  hence  q  =  2  in  (2)).  A  DSOCP  problem  in 
primal  standard  form  (see  [1])  is 


and  its  dual  problem  (see  [1]) 


min  cTx 

s.t.  Ax  =  b 

x  hr  0 

max  bJ  y 

s.t.  ATy  +  z  =  c 

z  hr  0 


(3) 


(4) 


where  A  G  Mmxn,  b  G  and  c  G  Rn  constitute  given  data,  x  G  Rn  is  the  primal  variable,  and 
y  G  Mm  and  z  G  Mn  are  the  dual  variables. 

Since  DSOCP  is  a  special  case  of  DMOCP,  all  problems  that  can  be  formulated  as  DSOCPs,  such 
as  DLPs  (the  DSOCP  problems  (3)  and  (4)  reduce  to  DLP  problems  when  r  =  n),  strictly  convex 
deterministic  quadratic  programs  (DQPs),  convex  quadratically  constrained  quadratic  programs 
(QCQPs),  and  problems  with  hyperbolic  constraints  (see  [1,  9])  are  special  cases  of  DSOCPs  and  of 
DMOCPs.  The  survey  paper  of  Lobo,  et  al.  [9]  discusses  DSOCPs  with  a  number  of  applications 
in  many  areas  including  a  variety  of  engineering  applications. 


2.1.1  Examples:  Norm  minimization  problems 

In  [1]  Alizadeh  and  Goldfarb  presented  DSOCP  formulations  of  three  norm  minimization  problems 
where  the  norm  is  the  Euclidean  norm.  In  this  subsection  we  show  how  extensions  of  these  three 
problems  where  we  use  arbitrary  p  norms  lead  to  DMOCPs.  Let  v%  =  AiX  +  G  Mni_1,  i  = 
1, 2,  •  •  •  ,  r.  The  following  norm  minimization  problems  can  be  cast  as  DMOCPs: 

1.  Minimization  of  the  sum  of  norms: 

The  problem  min  Ya= i  ll^llpi  can  be  formulated  as 

min  ELi  U 

s.t.  AiX  +  bi  =  Vi,  z  =  1,2,---  ,  r 

(ti;vi;t2-,v2;---  ]tr',vr )  —(pltp2r--,pr)  ^ 
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2.  Minimization  of  the  maximum  of  norms: 

The  problem  min  maxi<,<r  ||u,;||Pt  can  be  expressed  as  the  DMOCP  problem 

min  t 

s.t.  AiX  +  bj  =  Vi,  i  =  1,  2,  •  •  •  ,  r 

; t-, vr)  y{pi^...  pr)  o 

3.  Minimization  of  the  sum  of  the  k  largest  norms: 

More  generally,  the  problem  of  minimizing  the  sum  of  the  k  largest  norms  can  also  be  cast  as 
DMOCPs.  Let  the  norms  ||v[i]||P[1],  ||v[2]llp[2],  ■ ' '  ,  ||v[r]||pM  be  the  norms  ||ui||Pl,  ||v2||p2,  •  •  •  ,  IKHPr 
sorted  in  nonincreasing  order.  Then  the  problem  Yli=i  ||'yp]|Li,  can  be  formulated  as  (see 
also  [1]  or  [9]  and  the  related  references  contained  therein) 

min  Y^i=i  si  +  ht 

s.t.  AiX  +  bi  =Vi,  z  =  1,2,---  ,  r 

(si  +  t-,v1;s2  +  t;v2',---  ;sr  +  t,vr)  h{pi)P2)... >Pr)  0 

Si  >  0 

2.2  Duality 

Since  DMOCPs  are  a  class  of  convex  optimization  problems,  we  can  develop  a  duality  theory  for 
them.  A  forthcoming  paper  [2]  presents  such  a  duality  theory.  Here  we  indicate  weak  and  strong 
duality  for  the  pair  (P,  D)  as  justification  for  referring  to  them  as  a  primal  dual  pair. 

We  first  show  that  the  dual  of  the  pth- order  cone  of  dimension  n  is  the  qth- order  cone  of 
dimension  n,  where  q  is  the  conjugate  to  p.  For  any  cone  /C,  the  dual  cone  K,*  is  defined  by 
1C*  :=  {y  G  Mn  :  yTx  >  0,\/x  G  /C}. 

Lemma  1  Qp*  =  Qq,  where  1  <  p  <  oo  and  q  is  the  conjugate  top.  More  generally,  Q(Pl,P2)...  ,Pr)*  = 
Q(qi,q 2,...  ,qr),  where  1  <  Pi  <  oo  and  qt  is  the  conjugate  to  pi  for  i  =  1, 2,  •  •  •  ,  r. 

Proof.  We  assume  that  Qp  C  K".  The  proof  of  the  second  part  trivially  follows  from  the  first 
part.  To  prove  the  first  part,  we  first  prove  that  Qq  C  Qp*.  Let  x  =  ( xq\x )  G  Qq,  we  show 
that  x  G  Qp*  by  verifying  that  xJy  >  0  for  any  y  G  Qp.  So  let  y  =  (yo',y)  G  Qp.  Then 
x Jy  =  xo  2/o  +  xTy  >  ||®||q||y||p  +  xJy  >  \xTy  \  +  x^y  >  0,  where  the  first  inequality  follows 
from  the  fact  that  x  G  Qq  and  y  G  Qp  and  the  second  one  from  Holder’s  inequality.  Now  we 
show  Qp*  C  Qq.  Let  y  =  (yo',y)  €  Qp*,  we  show  that  y  G  Qq  by  verifying  that  yo  >  \\y\\q.  This 
is  trivial  if  y  =  0  or  p  =  oo.  If  y  ^  0  and  1  <  p  <  oo,  let  u  :=  (yip^q',y2p^q',  •  ■  •  ;  yn-\p^q)  and 
consider  x  :=  (||u||p;  —u)  G  Qp.  Then  by  using  Holder’s  inequality,  where  the  equality  is  attained, 
we  obtain  0  <  xTy  =  ||tt||p  yo  ~  uJy  =  ||w||p  yo  —  IMIP||i/||<j  =  ||u||p  (yo  —  ||y||9).  This  gives  that 

yo>\\y\\q-  n 

It  follows  from  this  lemma  that  the  second-order  cone  is  self-dual,  i.e.,  Q2*  =  Q2-  From  this 
lemma  we  also  deduce  that  the  p^-order  cone  is  reflexive,  i.e.,  Qp**  =  Qp,  and  more  generally,  also 
Q(pi,P2:.”  :pr)  is  reflexive.  On  the  basis  of  this  fact,  it  is  natural  to  infer  that  the  dual  of  the  dual  is 
the  primal.  Using  the  above  lemma,  we  can  prove  the  following  weak  duality  property. 
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Theorem  1  (Weak  duality)  If  x  is  any  primal  feasible  solution  of  (P)  and  ( y ,  z)  is  any  dual  feasible 
solution  of  (D),  then  the  duality  gap  cTx  —  bT y  =  xTz  >  0. 

Proof.  Note  that  cTx—bTy  =  (AJy+z)Jx  —  bJy  =  yT Ax+zJx  —  yJb  =  yJ(Ax  —  b)+zTx  =  xT z. 
Since  x  G  Q(Pl.P2....  ,Pr)  and  z  G  Q{qi,q2,  ~  ,qr)  =  Q{pi,p2,-  ,Pr)*>  we  conclude  that  xT z  >0.  □ 

We  can  now  apply  the  duality  relations  [11,  Theorem  4.2.1]  to  obtain: 

Theorem  2  (Strong  duality)  Assume  the  interior  of  (P)  to  be  not  empty  and  the  objective  value 
of  (P)  to  be  bounded  below  on  the  feasible  region.  Then  (D)  is  solvable  and  the  optimal  objective 
value  p*  of  (P)  and  the  dual  objective  value  d*  satisfy  the  relation  p*  =  d*. 

3  Definition  of  an  SMOCP 

In  this  section  we  define  two-stage  stochastic  multi-order  cone  programs  (SMCOPs)  with  recourse 
based  on  DMOCP  (P)  analogous  to  the  way  SLPs  is  defined  based  on  DLPs.  Let  ri,r2  >  1  be 
integers.  For  i  =  1,2,  ■  ■  ■  ,r\  and  j  =  1,2,  •••  ,r2,  let  Pu,P2j  £  [Loo]  and  mi,m2,n\,ri2,n\i,n2j 
be  positive  integers  such  that  ni  =  Y^i=inH  and  n-2  =  X][=i  n2j-  An  SMOCP  with  recourse  in 
primal  standard  form  is  defined  based  on  deterministic  data  A  G  b  £  Rmi  and  c  G  l"1  and 

random  data  T  G  Wn2'Kni,W  g  Rm2Xn2,h  G  Mm2  and  d  G  Mn2  whose  realizations  depend  on  an 
underlying  outcome  ui  in  an  event  space  0  with  a  known  probability  function  P.  Given  this  data, 
an  SMOCP  with  recourse  in  primal  standard  form  is 

min  cJx  +  E  [Q(x,  a;)] 

s.t.  Ax  =  b  (5) 

x>- ,  x0 

where  x  G  M"1  =  l”11  x  l"12  x  •  •  •  x  Mniri  is  the  first-stage  decision  variable  and  Q(x,co)  is  the 
minimum  of  the  problem 


min  d(co)Ty 

s.t.  T(w)x  +  W{oS)y  =  h(u)  (6) 

^  “  {P2 1  iP22  > " '  ,P2r2)^ 

where  y  G  Mn2  =  l”21  x  Rn22  x  •  •  •  x  R”2r2  is  the  second-stage  variable  and 

E [Q(x,u)\  :=  /  Q(x,Lo)P(dio). 

Jn 


3.1  Special  cases  of  SMOCPs 

In  this  part  we  present  some  important  special  cases  of  SMOCPs.  Stochastic  pth-order  cone  pro¬ 
grams  (SPOCPs)  are  a  special  case  of  SMOCPs  which  occurs  when  pu  =  p2j  =  p  >  1  for  all 
z  =  1,2,---  ,  r\  and  j  =  1,  2,  •  •  •  ,  r2  in  (5,  6).  An  SPOCP  problem  therefore  is 

min  cTx  +  E  [Q(x,  cu)] 

s.t.  Ax  =  b  (7) 
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where  x  G  Rni  =  R”11  x  Rni2  x  •  •  •  x  Rniri  is  the  first-stage  decision  variable  and  Q(x,oj)  is  the 
minimum  of  the  problem 


min  d(oj)Ty 

s.t.  T(w)x  +  W(uj)y  =  h(oj)  (8) 

y >-  ,  x0 
y  —r2{p) 

where  y  G  Rn2  =  R”21  x  Rn22  x  •  •  •  x  Rra2r2  is  the  second-stage  variable  and 

E[Q(a:,  cu)]  :=  f  Q(x,uj)P(dco). 
io 

Stochastic  second-order  cone  programs  (SSOCPs)  are  a  special  case  of  SPOCPs  (and  hence  of 
SMOCPs)  which  occurs  when  p  =  2  in  (7,  8).  An  SSOCP  problem  (see  also  [10])  is 

min  cTx  +  E  [Q(x,  cu)] 

s.t.  Ax  =  b  (9) 

xPri0 

where  x  G  Rni  =  Rnu  x  l"12  x  •  •  •  x  Rrairi  is  the  first-stage  decision  variable  and  Q(x,co)  is  the 
minimum  of  the  problem 


min  d(co)Ty 

s.t.  T(w)x  +  W(u)y  =  h(u)  (10) 

y^r2  o 

where  y  G  Mn2  =  M"21  x  Rn22  x  •  •  •  x  R”2r2  is  the  second-stage  variable  and 

E [Q(x,uj)]  :=  /  Q(x,co)P(dco). 

Jn 

Note  that  when  r\  =  ri\  and  r2  =  n2,  SSOCP  (9,  10)  reduces  to  an  SLP.  So  SLPs  are  a  special 
case  of  SSOCPs. 

In  the  rest  of  this  section,  we  will  show  that  stochastic  quadratic  programs  (SQPs)  can  also 
be  cast  as  SSOCPs.  Our  proof  is  parallel  to  the  proof  of  the  fact  that  DQPs  is  a  subclass  of 
DSOCPs  (see  [1]).  Recall  that  a  two-stage  SQP  (with  recourse)  is  defined  based  on  deterministic 
data  C  G  R"lV  ni,C  0,  c  G  Rni,7l  G  Mmixni  and  b  G  Rmi;  and  random  data  H  G  Rn2V  n2,H  y 
0 ,d  G  R”2,T  G  Rm2 xni .  W  G  Rm2Xn2,  and  h  G  R"12  whose  realizations  depend  on  an  underlying 
outcome  in  an  event  space  fl  with  a  known  probability  function  P.  Given  this  data,  an  SQP  with 
recourse  is 


min  qi(x,ui)  =  xTCx  +  cTx  +  E[Q(aj,  a;)] 

s.t.  Ax  =  b  (11) 

x  >  0 

where  x  G  Rni  is  the  first-stage  decision  variable,  Q(x,lo)  is  the  minimum  of  the  problem 

min  q2(y,uj)  =  yT  H(w)y  +  d(u)J  y 

s.t.  T(u)x  +  W{uj)y  =  h{ij)  (12) 

y  >  o 
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where  y  €  M™2  is  the  second-stage  variable,  and 


E [Q{x,u)\  :=  /  Q(x,u>)P(du>). 


Observe  that  the  objective  function  of  (11)  can  be  written  as  (see  §2  in  [1]), 
qi(xi,ui)  =  ||w||2  +  E [Q(x,u)\  —  -^cTC~1c  where  u  =  C^^x  + 

Similarly,  the  objective  function  of  (12)  can  be  written  as 

<72(2/, xi)  =  IlDll2  —  i  d(co)T H(uj)~ld  (u>)  where  v  =  {uj)-1  ^  2  d{u) . 


Thus,  problem  (11,  12)  can  be  transformed  into  the  SSOCP: 

min  no 

s.t.  u  —  C*1/2*  =  C-1!2* 

Ax  =  b 

x  >  0,  (no;  u)  0 

where  Q(x,uj)  is  the  minimum  of  the  problem 


(13) 


mm  no 

s.t.  v  —  H{uj)l/2y  =  ^  H (u))-1  / 2 d(ijj) 
T(uS)x  +  W(u>)y  =  h(  to) 
y  >  0;  (u0;v)  t  0 


(14) 


where 


E [Q(x,u)\  :=  /  Q{x,Lo)P(dio). 

Note  that  both  problems  (the  SQP  problem  and  the  SSOCP  problem)  will  have  the  same 
minimizers,  but  their  optimal  objective  values  are  equal  up  to  constants.  More  precisely,  the 
difference  between  the  optimal  objective  values  of  (12)  and  (14)  would  be  —  \d(u:)J PL(oj)~l  d{ u). 
Consequently,  the  optimal  objective  values  of  (11,  12)  and  (13,  14)  will  differ  by 


--  cTC,-1c-  - 
2  2 


(d(cn)T  if(cn)_1d(cn)^  P(dw). 


4  Definitions  of  a  DMIMOCP  and  a  O-IDMOCP 

In  this  section  we  introduce  two  important  related  problems  that  result  when  decision  variables 
in  an  MOCP  can  only  take  integer  values.  Consider  the  DMOCP  problem  (P).  If  we  require  an 
additional  constraint  that  a  subset  of  the  variables  have  to  attain  0-1  values,  then  we  are  interested 
in  optimization  problem  of  the  form 


min  c^x 
s.t.  Ax  =  b 

x>~ ,  ,0 

—  (Pl,P2,---,Pr) 

xk  G  {0, 1},  k  G  r 

where  T  C  {1,  2,  •  •  •  ,n},  the  decision  variable  x  G  R”  has  some  of  its  components  Xk  (k  G  T)  with 
integer  values  and  bounded  by  ak,Pk  £  R-  This  class  of  optimization  problems  may  be  termed  as 
deterministic  0-1  multi-order  cone  programs  (O-lDMOCPs). 

A  more  general  and  interesting  problem  when  in  an  DMOCP  some  variables  can  only  take 
integer  values.  If  we  are  given  the  same  data  A,  b,  and  c  as  in  (P),  then  we  are  interested  in  the 
problem  of  the  form 


min 

s.t. 


CTX 

Ax  =  b 


X>z 


(Pl>P2,-,Pr) 


0 


Xk  G  [ak,/3k]f]Z,  k  €  r 


(15) 


where  T  C  {1,  2,  •  •  •  ,  n},  the  decision  variable  x  G  R”  has  some  of  its  components  Xk  ( k  G  T)  with 
integer  values  and  bounded  by  otk,Pk  £  R-  This  class  of  optimization  problems  may  be  termed 
as  deterministic  mixed  integer  multi-order  cone  programs  (DMIMOCPs).  The  relationships  among 
DMIMOCPs,  deterministic  mixed  integer  pth-order  cone  programs  (DMIPOCPs)  (which  occurs 
when  pi  =  p  >  1  for  all  i  =  1,2,  ••  •  ,  r),  deterministic  mixed  integer  second-order  cone  programs 
(DMISOCPs)  [7]  (which  occurs  when  pi  =  2  for  all  i  =  1,  2,  •  ■  ■  ,r),  and  deterministic  mixed  integer 
linear  programs  (DMILCPs)  (or  deterministic  mixed  integer  quadratic  programs)  are  the  same  as 
those  among  DMOCPs,  DPOCPs,  DSOCPs,  and  DMILPs  (or  DMIQPs),  respectively. 

We  can  also  handle  uncertainty  in  data  defining  DMIMOCPs  by  defining  two-stage  stochastic 
mixed  integer  multi-order  cone  programs  (SMIMCOPs)  with  recourse  (which  generalizes  two-stage 
stochastic  mixed  integer  second-order  cone  programs  [3])  based  on  DMIMOCP  (15)  analogous  to 
the  way  SMOCPs  (5,  6)  are  defined  based  on  DMOCPs  (P).  See  Figure  1  which  shows  conceptual 
relationships  among  the  optimization  problems  over  multi-order  cones  described  above  and  their 
special  cases. 


5  An  application 

Our  application  is  four  versions  of  the  facility  location  problem  (FLP).  For  these  four  versions  we 
present  problem  descriptions  leading  to  a  DMOCP  model,  an  SMOCP  model,  a  O-IDMOCP  model, 
and  a  DMIMOCP  model. 

In  FLPs  we  are  interested  in  choosing  a  location  to  build  a  new  facility  or  locations  to  build 
multiple  new  facilities  so  that  an  appropriate  measure  of  distance  from  the  new  facilities  to  existing 
facilities  is  minimized.  FLPs  arise  when  decisions  on  locating  airports,  regional  campuses,  wireless 
communications  towers,  etc.  are  to  be  made.  There  are  different  ways  of  classifying  FLPs.  Following 
are  some  of  these  ways  (see  also  [14]): 

•  We  can  classify  FLPs  based  on  the  number  of  new  facilities  in  the  following  sense:  if  we 
add  only  one  new  facility  then  we  get  a  problem  known  as  a  single  facility  location  problem 
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Figure  1:  Conceptual  relationships  among  the  optimization  problems  over  multi-order  cones  and 
their  special  cases. 


(SFLP),  while  if  we  add  multiple  new  facilities  instead  of  adding  only  one,  then  we  get  a  more 
general  problem  known  as  a  multiple  facility  location  problem  (MFLP). 

•  Another  way  of  classification  is  based  on  the  distance  measure  used  in  the  model  between  the 
facilities.  If  we  use  the  Euclidean  distance  then  these  problems  are  called  Euclidean  facility 
location  problems  (EFLPs),  if  we  use  the  rectilinear  distance  (also  known  as  L\  distance, 
city  block  distance,  or  Manhattan  distance)  then  these  problems  are  called  rectilinear  facility 
location  problems  (RFLPs).  Furthermore,  in  some  applications  we  use  both  the  Euclidean 
and  the  rectilinear  distances  (based  on  the  relationships  between  the  pairs  of  facilities)  as 
the  distance  measures  used  in  the  model  between  the  facilities  to  get  a  mixed  of  EFLPs  and 
RFLPs  that  we  refer  to  as  Euclidean-rectilinear  facility  location  problems  (ERFLPs). 

•  When  the  new  facilities  can  be  placed  any  place  in  solution  space,  the  problem  is  called  a 
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continuous  facility  location  problem  (CFLP),  but  usually  the  decision  maker  needs  the  new 
facilities  to  be  placed  at  specific  locations  (called  nodes)  and  not  in  any  place  in  the  solution 
space.  In  this  case  the  problem  is  called  a  discrete  facility  location  problem  (DFLP). 

•  In  some  applications,  the  locations  of  existing  facilities  cannot  be  fully  specified  because  the 
locations  of  some  of  them  depend  on  information  not  available  at  the  time  when  decision  needs 
to  be  made  but  will  only  be  available  at  a  later  point  in  time.  In  this  case,  we  are  interested  in 
stochastic  facility  location  problems  (or  abbreviated  as  stochastic  FLPs).  When  the  locations 
of  all  old  facilities  are  fully  specified,  FLPs  are  called  deterministic  facility  location  problems 
(or  abbreviated  as  deterministic  FLPs). 

FLPs  have  seen  a  great  deal  of  recent  research  activity.  For  further  details,  consult  the  book 
of  Tompkins  and  et  al.  [14].  In  particular,  deterministic  Euclidean  facility  location  problems  are 
often  cited  as  an  application  of  deterministic  second-order  cone  programs  (see  for  example  [15]  and 
[8]).  Each  one  of  the  next  subsections  is  devoted  to  a  version  of  ERFLPs.  Specifically,  we  consider 
deterministic  continuous  Euclidean-rectilinear  facility  location  problems  (deterministic  CERFLPs) 
which  leads  to  a  DMOCP  model,  stochastic  continuous  Euclidean-rectilinear  facility  location  prob¬ 
lems  (stochastic  CERFLPs)  which  leads  to  an  SMOCP  model,  deterministic  discrete  Euclidean- 
rectilinear  facility  location  problems  (deterministic  DERFLPs)  which  leads  to  a  O-IDMOCP  model, 
and  deterministic  ERFLPs  with  integrality  constraints  which  leads  to  a  DMIMOCP  model. 

5.1  Deterministic  CERFLPs — A  DMOCP  model 

In  deterministic  single  ERFLPs,  we  are  interested  in  choosing  a  location  to  build  a  new  facility 
among  existing  facilities  so  that  this  location  minimizes  the  sum  of  a  weighted  (either  Euclidean  or 
rectilinear)  distance  to  all  existing  facilities. 

Assume  that  we  are  given  r  +  s  existing  facilities  represented  by  the  fixed  points  oi,  <Z2,  •  •  •  ,  ar, 
ar+ 1,  ar+ 2,  ■  •  • o,r+s  in  Rn,  and  we  plan  to  place  a  new  facility  represented  by  x  so  that  we  minimize 
the  weighted  sum  of  the  Euclidean  distances  between  x  and  each  of  the  points  Oi,  02,  •  •  ■  ,  ar  and 
the  weighted  sum  of  the  rectilinear  distances  between  x  and  each  of  the  points  ar+ 1,  ar+ 2,  •  •  •  ,  ar+s. 
This  leads  us  to  the  problem 

min  Yh=i  wi  \\x  ~  1 1 2  +  Ya=?+  iwi  \\x  ~  ®i||i 

or,  alternatively,  to  the  problem 

min  X)i=i  wi  U 

s.t.  d»i j  ^  ^  x  0 

(tr_l_i5  x  ar- |_i5  5  tr-\-s'i  'E  ^r+s)  ^»<i>  0 

where  Wi  is  the  weight  associated  with  the  zth  existing  facility  and  the  new  facility  for  i  = 
1,2, ...  ,r  +  s. 

In  deterministic  multiple  ERFLPs  we  add  m  new  facilities,  namely  x\,  *2,  •  •  •  ,  xm  G  R",  instead 
of  adding  only  one.  We  have  two  cases  depending  whether  or  not  there  is  an  interaction  among 
the  new  facilities  in  the  underlying  model.  If  there  is  no  interaction  between  the  new  facilities,  we 
are  just  concerned  in  minimizing  the  weighted  sums  of  the  distance  between  each  one  of  the  new 
facilities  and  each  one  of  the  fixed  facilities.  In  other  words,  we  solve  the  following  DMOCP  model: 
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Em  sr~^r-\-s  « 

j= 1  Z)t=l  wij  tij 

S.t.  (ilj  5  'Ej  5***5  trj  5  ^r*)  ^Zr(2)  0,  J  —  1,  2,  •  •  •  ,  7TI 

(*(r+l)j>  ^r+1?  *  *  *  ?  ^(r+s)j>  ^r+s)  ^s(l)  J  1,  2,  •  •  •  ,  777/ 


(16) 


where  w^j  is  the  weight  associated  with  the  zth  existing  facility  and  the  jth  new  facility  for  j  = 
1,2, ...  ,m  and  i  =  1, 2, . . . ,  r  +  s. 

If  interaction  exists  among  the  new  facilities,  then,  in  addition  to  the  above  requirements,  we 
need  to  minimize  the  sum  of  the  (either  Euclidean  or  rectilinear)  distances  between  each  pair  of 
the  new  facilities.  Let  1  <  l  <  m  and  assume  that  we  are  required  to  minimize  the  weighted  sum 
of  the  Euclidean  distances  between  each  pair  of  the  new  facilities  *i,  *2,  •  •  •  ixi  and  the  weighted 
sum  of  the  rectilinear  distances  between  each  pair  of  the  new  facilities  cc;_|_i,  cc;+2>  ■  •  •  ,  xm.  In  this 
case,  we  are  interested  in  a  model  of  the  form: 


min 

s.t. 


Em  v-vr+s 

j= 1 


Wij 


{t\j  \  Xj  d\ 5  •  • 


tij  +  2  Xff =1  Wjj'tjj' 

i  trji  xj  ®r)  (Hr('2)  0,  j  =  1,  2,  •  •  • 


,  m 


(^(r+l)j>  —  ar+li 

wjO’+i)’  —  ®j+ii 


i  ^(r+s)ji  Xj  —  ttr+s)  — s(l)  7  1,  2,  •  •  ■  ,  771 

;  tji ;  —  a?/)  ^I(j— j)(2)  0,  j  =  1, 2,  •  •  •  ,  l  —  1 

j  tjm\  Xj  —  Xm )  ^(m— j)(l)  ^5  J  =  ^  +  I5  2,  •  •  • 


m  —  1 


(17) 


where  Wjj>  is  the  weight  associated  with  the  new  facilities  j1  and  j  for  j'  =  1,2, ...  ,j  —  1  and 

j  =  2, 3, . . . ,  m. 

5.2  Stochastic  CERFLPs — An  SMOCP  model 

Before  we  describe  the  stochastic  version  of  this  generic  application,  we  indicate  a  more  concrete 
version  of  it.  Assume  that  we  have  a  new  city  with  many  suburbs  and  we  want  to  build  a  hospital 
for  treating  the  residents  of  this  city.  Some  people  live  in  the  city  at  the  present  time.  As  the 
city  expands,  many  houses  in  new  suburbs  need  to  be  built  and  the  locations  of  these  suburbs  will 
be  known  in  the  future.  Our  goal  is  to  find  the  best  location  of  this  hospital  so  that  it  can  serve 
the  current  suburbs  and  the  new  ones.  This  location  must  be  determined  at  the  current  time  and 
before  information  about  the  locations  of  the  new  suburbs  become  available. 

Generally  speaking,  let  a,2,  ■  ■  ■  ,  ari,  ari+i,  ari+2,  •  •  •  ,  ari+Sl  be  fixed  points  in  W1  represent¬ 
ing  the  coordinates  of  n+si  existing  fixed  facilities  and  di(w),  a 2(0;),  •  •  •  ,  ar2( u>),  dr2+i(u;),  dr2+2(w) 
,  •  •  •  ,  ar2+S2  (c j)  be  random  points  in  Mn  representing  the  coordinates  of  +  s 2  random  facilities 
whose  realizations  depends  on  an  underlying  outcome  u>  in  an  event  space  P  with  a  known  proba¬ 
bility  function  P. 

Suppose  that  at  present  we  do  not  know  the  realizations  of  r2  +  s 2  random  facilities,  and  that 
at  some  point  in  time  in  future  the  realizations  of  these  r2  +  s 2  random  facilities  become  known. 
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Our  goal  is  to  locate  m  new  facilities  X\,  X2,  -  -  *  ,  xm  G  that  minimize  the  the  following  sums: 

•  The  weighted  sums  of  the  Euclidean  distance  between  each  one  of  the  new  facilities  and  each 
one  of  the  fixed  facilities  oi,  <22,  •  •  •  ,  ari, 

•  The  weighted  sums  of  the  rectilinear  distance  between  each  one  of  the  new  facilities  and  each 
one  of  the  fixed  facilities  ari+i,  ari+2,  •  •  •  ,  ari+si, 

•  The  expected  weighted  sums  of  the  Euclidean  distance  between  each  one  of  the  new  facilities 
and  each  one  of  the  random  facilities  di(u;),  <22 [oS)t  ■  ■  ■  ,  ar2(u ;), 

•  The  expected  weighted  sums  of  the  rectilinear  distance  between  each  one  of  the  new  facilities 
and  each  one  of  the  random  facilities  dr2+i(u;),  dr2+2(u;)  ,  •  •  •  ,  ar2+S2( uo). 

Note  that  this  decision  needs  to  be  made  before  the  realizations  of  the  r2  +  S2  random  facilities 
become  available.  If  there  is  no  interaction  between  the  new  facilities,  then  the  DMOCP  model 
(16)  becomes  the  following  SMOCP  model: 

min  YJjL  1  Wij  **3  +  E  [Q(x  1;  *  *  *  ;  xm,u)} 

S.t.  ( t\j  5  Xj  CL\ 5  *  *  *  5 1 rij  5  xj  ®ri )  (2)  3  =  1?  2,  *  *  *  ,  777/ 

(t(ri+l)j  ’  ®ri+l  j  *  *  *  ?  ^(ri+si)j  5  ^ri+si )  (1)  0,  J  1,  2,  •  •  •  ,  777- 

where  Q(a?i;  •  •  •  ;  xm,cj)  is  the  minimum  of  the  problem 

min  E^lr  M  Uj 

s.t.  -  «iM;  •  •  •  ;tr2j;Xj  -  ar2(u ))  ^r2<2)  0,  j  =  1,2,  •••  ,m 

(tlra+llj,  ®(r2+l)  (^)j  >  ^(r2+S2)j)  ®(r2+S2)  (^))  —  S2(l)  7  '  '  '  i™ 


1E[Q(*1,  •  •  •  j®mi  k-0]  • —  /  i  Xm,  Uj')P(duj) . 

Jci 


where  Wij  is  the  weight  associated  with  the  ith  existing  facility  and  the  jth  new  facility  for  j  = 
1,2, ...  ,7ii  and  i  =  1,2, ... ,  n  +  si,  and  Wij(  uj)  is  the  weight  associated  with  the  ith  random  existing 
facility  and  the  jth  new  facility  for  j  =  1,2, ...  ,m  and  i  =  1,  2, . . . ,  72  +  s 2. 

If  interaction  exists  among  the  new  facilities,  then  the  DMOCP  model  (17)  becomes  the  following 
SMOCP  model: 
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min 

s.t. 


Ej=i  E[itsi  wv  +  E™=2  £7=\  Wjj'tjf  +  E  [Q(*i;  •  ■  ■ ;  *m, ^)] 

(tij  5  Q-i  5***5  trij  5  ®ri  )  (2)  ®5  J  =  1?  2,  •  •  •  ,  Vfl 


(*(n+l)j  ’  ®ri+l  j  *  *  *  5  ^(n+si)j  ’  ^r*i+si )  ^si (1)  J  1?  2,  *  *  *  ,  Tfl 

V'jij+i)  ’  ®?+i»  5  ^ jh  *^j  *^/)  —  (/— i)( 2)  .7  1?  2,  •  •  •  ,  l  1 

(^j'(j+l)  5  ®7+l?  5  tjm]  Xm  )  h{m-j){ i)  0,  j  =  l  +  1, 2,  ■  •  •  ,  m  -  1 


where  Q(x i;  •  •  •  ;  xm,ui )  is  the  minimum  of  the  problem 


min 

s.t. 


£7=i  £SS2  ^  M  4-  +  £7=2  £7=i  £/0/ 

-Oi(o;);--  -  ;tr2j;Xj  -  ar2(u>))  hr2(2)  0,  j  =  1,2,  •••  ,m 

(t(r2+l)jiXj  ~ '  ®(r2+l)  (w)>  i  ^(r2+s2)j  i  Xj  ~  a{r2+s2)(.1^'))  CZs2(i)  0)  J  1)2,  •••  ,777. 

(Ai(i+i)’  ®j+i!  )  tjh  xj  xi)  — (i— j)(2)  0)  j  1)  2)  ■  *  *  ,  l  1 

(7/(i+ 1)>  ®j+ij  !  ijm]  Xj  xm  )  ^(m-j)(i>  0,  j  =  l  +  1, 2,  •  •  •  ,  m  -  1 


and 


E[1?(*®1 )  ‘  • —  /  Q(®lj  iXm,0j)P{dlo). 

Jn 


where  utjji  is  the  weight  associated  with  the  new  facilities  j'  and  j  for  j 1  =  1,2, ...  ,j  —  1  and 

J  =  2, 3, . . . ,  m. 


5.3  Deterministic  DERFLPs — A  O-IDMOCP  model 

We  consider  the  discrete  version  of  the  problem  by  assuming  that  the  new  facilities  aq,  aq,  •  **  , 
need  to  be  placed  at  specific  locations  and  not  in  any  place  in  2-  or  3-  (or  higher)  dimensional 
space.  Let  the  points  V\,V2,  •  •  •  ,  vk  €  Mn  represent  these  specific  locations  where  k  >  m.  So,  we 
add  the  constraint  xt  e  {v\,V2,  •  •  •  ,  r’fc}  for  i  =  1, 2,  ■  ■  •  ,  m.  Clearly,  for  z  =  1,2,---  ,  m,  the  above 
constraint  can  be  replaced  by  the  following  linear  and  binary  constraints: 


Xi  =  vx  ya  +v2yi2~\ - \~vk  yik, 

Vi l  +Vi2-\ - t  Vik  =  1,  and 

Vi  =  (m;  Vii\ •  •  •  ;  y%k)  e  {o,  i}fc. 

We  also  assume  that  we  cannot  place  more  than  one  facility  at  each  location.  Consequently,  we 
add  the  following  constraints: 

(i;  2/iz;  2/2/;  -  -  -  vymi)  £i)  o,  for  i  =  1, 2,  —  ,  k. 
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If  there  is  no  interaction  between  the  new  facilities,  then  the  DMOCP  model  (16)  becomes  the 
following  O-IDMOCP  model: 

Era  # 

j  ;  £*=i  wij  Uj 

s.t.  (ti jj  Xj  cl\ ,  ,  trj ,  Xj  ctr)  >rr(2)  0,  j  —  1,  2,  •  •  •  ,  m 

(Z(r+l)jp  xj  ®r+l!  !  t(r+s)ji  xj  ®r+s)  ^Zs(l)  0,  J  1)  2,  ■  •  •  771 

=  «i  2/ii  +  v2  Vi2  H - b  Vik,  *  =  1,2,  —  ,  m 

(i; 2/iz; •  •  •  ;ymi)  t<i>  o,  for  z  =  1, 2, - ••  ,k 
1T  t/j  =  1,  Hi  G  {0,  l}fe,  *  =  1,2,---  ,  m 

If  interaction  exists  among  the  new  facilities,  then  the  DMOCP  model  (17)  becomes  the  following 
O-IDMOCP  model: 


min 

s.t. 


Era  v^r+s 

j  =  l  Li=l  Wij 

(t\j  Q-i  5  * 


tij  +  Ylj=2  Wjj'tjj' 

i  trj  i  Xj  flr)  — r(2)  0,  J  =  1,  2,  •  •  • 


,  m 


(t(r_^_  i  jj ,  Xj  ,  Z(r+s)j)  xj  ^r+s)  >IS<1>  Hi  J  1,2,"""  ,  TYl' 

(*j(j+l)i  —  *7+1)  i  tjli  xj  ~  xl )  — (J— j)( 2)  H,  j  =  1,  2,  •  •  •  ,  l  —  1 

(^i(i+l)’  —  *j+l>  j  tjm'i  xj  —  xm)  — (m— j)(l)  0,  3  =  Z  +  1,  2,  •  •  •  ,  777.  —  1 

=  Vi  y*i  +  ^2  Vi2  H - 1-  Vk  yik,  *  =  1,2,---  ,  m 

(i;  2/iz; */20  •  -  -  ;ymj)  ^<i)  o,  for  /  1.2.  •••,./,• 

1T  y*  =  1,  Vi  €  {0,  l}fc,  *  =  1,2,---  ,  m 


For  l  =  1,2,  ■■■  ,  k,  let  27  =  1  if  the  location  V{  is  chosen,  and  0  otherwise.  Then,  we  can  go 
further,  and  consider  more  assumptions:  Let  ki,k2,  ^'3,  ^'4  6  [1,/c]  be  integers  such  that  £+  <  k2 
and  &3  <  ZC4.  If  we  must  choose  at  most  &+  of  the  locations  Vi,V2,'"  , Vfc2,  then  we  impose  the 
constraints: 

{k\\Zi]Z2V  ■  ■  ;zk 2)  ^(1)  0,  and  2  G  {0,  l}fc. 

If  we  must  choose  at  most  k\  of  the  locations  vi,v2,---  , vk2 ,  or  at  most  k%  of  the  locations 
V\,  v2,  ■  ■  ■  ,  vk4,  then  we  impose  the  constraints: 

(kif;zi;z2;---  ;zk2)  0,  (fc3(l  -  f);  zi;  z2;  ■  ■  ■  ;zki)  0,  2  G  {0,  l}fc,  and  /  G  {0,1}. 
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5.4  Deterministic  ERFLPs  with  integrality  constraints — A  DMIMOCP  model 

In  some  problems  we  may  need  the  locations  to  have  integer-valued  coordinates.  In  most  cities, 
streets  are  laid  out  on  a  grid,  so  that  city  is  subdivided  into  small  numbered  blocks  that  are  square 
or  rectangular.  In  this  case,  usually  the  decision  maker  needs  the  new  facility  to  be  placed  at  the 
corners  of  the  city  blocks.  Thus,  for  each  i  6  Ac  {1,  2,  •  •  •  ,  m},  let  us  assume  that  the  variable  xt 
lies  in  the  hyper-rectangle  3”  =  {*,  :  Ci  <  Xi  <  77*,  Ci  £  Mn,  £  Kn}  and  has  to  be  integer-valued, 
i.e.  xt  must  be  in  the  grid  E”  P|  Zn.  Thus,  if  there  is  no  interaction  between  the  new  facilities,  then 
instead  of  solving  the  DMOCP  model  (16),  we  solve  the  following  DMIMOCP  model: 


min 


Em  x~^r-\-s  I 

7=1  J2i=  1  wij  Uj 


+  Er=3Ejr=>«4-' 


S.t.  (tl ji  Xj  d\,  ,  trj ,  ®r)  ^Zr( 2)  0,  J  —  1>  2,  ■  ■  ■  ,  771 

(^(r+l)j!  ~  °,T'+lj  j  ~  ar+s)  ^s(l)  J  =  I5  2,  •  •  •  ,  771 


£HgnZn,  fc£  A 


If  interaction  exists  among  the  new  facilities,  then  instead  of  solving  the  DMOCP  model  (17), 
we  solve  the  following  DMIMOCP  model: 


min 

s.t. 


7=1  Li=l  wij  tij  +  l^j=2  Ljfcl 


(t\j  ,  ®1,  j  trj ,  aij  Q.r)  ^j-^2)  0)  J  —  1)  2,  •  •  •  ,  771 

(^(r+l)ji  —  ®r+lj  >  t(r+s)ji  ®r+s)  CZs(l)  O5  7  I5  2,  ,  777. 

(^i(i+i)’  ®j+i)  !  tji ,  Xj  ai/)  ^(i— j)(2)  J  U  2,  •  •  •  5 /  1 

(*j(j+l)i  *1  —  *5+1)  i  tjm]  Xj  —  Xm)  ^(m— 0)  3 =  ^  T  1,  2,  •  •  •  ,  771  —  1 


£  32  nzB,  *  e  A 


Finally,  we  mention  that  it  is  useful  to  consider  the  stochastic  version  of  DERFLPs  to  obtain  a  0- 
1SMOCP  model,  and  to  consider  also,  parallelly,  the  stochastic  version  of  ERFLPs  with  integrality 
constraints  to  obtain  an  SMIMOCP  model. 

6  Concluding  remarks 

We  conclude  the  paper  by  indicating  possible  directions  for  future  research. 

As  indicated  in  Figure  1,  the  new  classes  of  optimization  problems  DPOCPs,  SPOCPs,  DMIPO- 
CPSs,  DMOCPs,  SMOCPs  and  DMIMOCPs  introduced  in  this  paper  are  natural  extensions  of 
deterministic,  stochastic  and  mixed  integer  second-order  cone  programs.  It  is  interesting  to  note 
that  semidefinite  programs  and  multi-order  cone  programs  both  include  second-order  cone  programs 


16 


as  a  special  case.  In  §5  we  presented  an  application  leading  to  multi-order  cone  programs.  It  is 
interesting  to  investigate  other  applicational  settings  leading  to  (deterministic,  stochastic  and  mixed 
integer)  multi-order  cone  programs.  Development  of  algorithms  for  such  multi-order  cone  programs 
which  in  turn  will  benefit  from  a  duality  theory  is  equally  interesting  and  important.  The  authors 
are  currently  exploring  these  research  directions. 
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